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Abstract. We present a generalization of the double semion topological quan¬ 
tum field theory to higher dimensions, as a theory of d — 1 dimensional surfaces 
in a d dimensional ambient space. We construct a local Hamiltonian which 
is a sum of commuting projectors and analyze the excitations and the ground 
state degeneracy. Defining a consistent set of local rules requires the sign struc¬ 
ture of the ground state wavefunction to depend not just on the number of 
disconnected surfaces, but also upon their higher Betti numbers through the 
semicharacteristic. For odd d the theory is related to the toric code by a local 
unitary transformation, but for even d the dimension of the space of zero energy 
ground states is in general different from the toric code and for even d > 2 it is 
also in general different from that of the twisted Z 2 Dijkgraaf-Witten model. 


1 . 

The notion of a topological quantum field theory (TQFT) had its initial success 
as a 2 -|- 1 dimensional theory employed by Witten to give the Jones polynomial a 
coordinate free interpretation. There are a few regimes in which TQFTs generalize 
straightforwardly to higher dimension, d + 1, d > 2. These include Dijkgraaf- 
Witten theories [1] or, more generally, theories based on the category of maps to a 
space Y with only hnitely many nontrivial homotopy groups which are, themselves, 
hnite [H]. Another family of examples are generalized toric codes (GTCs): Given 
0 < p < d there is a d-|-1 dimensional TQFT whose Hilbert space on a closed (not 
necessarily oriented) d-manifold X is spanned by the elements of Hp (A; Z 2 ). See 
Ref. [2] for the case p = 1, d = 2, see Ref. [3] for the case p = 1, d = 3, and see 
Ref. |1] for the case p = 2, d = 4. 

In dimension 2 -|- 1 the toric code is often studied hand-in-hand with a sister 
theory “double semions” or DS. In this theory the ground state vectors are multi¬ 
loops, with a loop fugacity of S = —1. The purpose of this paper is to explain 
(only) in the case p = d — 1 how the DS theory can be generalized to d > 2. We 
call this the generalized (double) semion theory (GDS); we parenthesize the word 
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“double” because these theories are not doubles in any sense that we know. We 
address: 

(1) gapped local Hamiltonian. This Hamiltonian will be a sum of commuting 
local terms, each of which is a projector. 

(2) The form of ground state wave functions—answering the question: What 
is the generalization of loop fugacity <5 = — 1? 

(3) The nature of loop and dual loop operations, which we now call “balloon” 
and “dual loop.’, and also the nature of excitations. 

(4) To what class of d-manifolds {X'^} should we restrict, to arrive at a theory 
still analogous to the GTCs, insofar as ground space degeneracies and 
balloon and loop operators are concerned. 

The multi-loops of the theory will be dehned on a hxed, hnite cellulation of 
some ambient manifold. Thus, the Hilbert space of the theory will be hnite di¬ 
mensional. Hamiltonians where the terms are commuting projectors are often 
studied as lattice model realizations of a TQFT. Such lattice model realizations 
are known in d = 2 via the Levin-Wen construction [5] and in d = 3 via the 
Walker-Wang construction [7]. Lattice models where the terms in the Hamiltonian 
commute are mathematically interesting as these properties simplify the mathe¬ 
matical treatment of these theories; further, stabilitiy of these Hamiltonians under 
small perturbations has been proven under fairly general assumptions|8l|9], so that 
results about topological order in these idealized Hamiltonians can often be shown 
to hold for more general Hamiltonians. 


2. Brief review of toric code (TC) and double semions (DS) 

The generalized toric code model takes as input a cellulation C of a d dimensional 
manifold X. There is a qubit degree of freedom for each p cell of C, for some given 
p which is also input to the theory, with the Hilbert space H being the tensor 
product of these Hilbert spaces (generalizations exists to qudit degrees of freedom; 
we do not consider this here). The Hamiltonian acting on H is a sum of two terms, 

(2.1) H = H+ + Ha, 
where 

( 2 . 2 ) H^= Ha= 

p — 1 cells e p + 1 cells c 

We dehne 


He = 



(2.3) 
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where Zf is the Pauli-Z operator on the cell / and df denotes the boundary of /. 
We dehne 

^2 ^ ~ n(d - 1) cells f &dc 

where the operator Wj is the Pauli-X operator acting on the cell /, changing the 
state of the cell. When we describe product of several Pauli-X operators on some 
set of cells later, we will sometimes refer to this as NOT on that set. 

These terms are all projectors, and they pairwise commute. The ground state 
degeneracy is given by ^ basis for the ground state wavefunctions of 

this theory is to take equal amplitude superpositions of closed p-chains in a given 
homology class. 

The theory is a topological quantum field theory for 0 < p < d. We now restrict 
to d = 2,p = 1. On a sphere, the ground state wavefunction of the toric code is a 
sum over equal amplitude superpositions of closed 1—chains. The double semion 
model on a sphere has a similar ground state wavefunction up to a sign factor: 
the ground state wavefunction is a sum over closed 1-chains, with an amplitude 
equal to —1 raised to a power equal to the number of closed loops. To define this 
amplitude unambiguously, one typically works on a hexagonal lattice (or other 
lattice with degree at most 3). If we instead consider a theory of loops on a square 
lattice, a “£gure-8” loop conhguration is possible for which it is ambiguous whether 
we have a single loop which crosses itself or two loops which touch at a corner. If 
the degree is at most 3 then closed 1-chains unambiguously dehne a disconnected 
set of closed loops: a spin j' on a 1-cell corresponds to loop being present on that 
1-cell and a spin j, corresponds to loop being absent. 

There are two distinct choices typically made for the double semion model 
Hamiltonian on a hexagonal lattice [S]. Both have the same ground state subspace 
and both have an excitation gap. In both, we again decompose H = + Hq, 

where H+ is the same as the toric code for p = 1, d = 2 and Hn = ^^2 cells c 
for He dehned as follows. In the hrst choice, one sets 


(2.5) 


H,= 




1 cells f G dc 


Xf 


where nuryic) 

the t state. Thus, H^dryic) jg equal to the product of the operators 


around 


is equal to the number of legs leaving the hexagon c which are in 

'1 0 ^ 

^0 

the legs leaving the hexagon. With this choice, the terms He commute with the 
terms He and the terms He pairwise commute with each other when restricted to 
the eigenspace of with vanishing eigenvalue. 

The second choice is to project the term in Eq. (2.5) into the zero eigenspace of 
the operators He for all vertices e in the hexagon c. With this choice, the terms 
all commute with each other. 
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Consider a configuration of loops on the sphere. The term He can induce a 
transition between two conhgurations differing by a surgery move: 

X^)( 

To be in the zero eigenspace of H, the two conhgurations related by this move 
must enter the wavefunction with the opposite sign. Given any conhguration of 
closed loops on the sphere, indeed this transition will always change the number 
of loops by ±1 and hence the sign rule for this transition is consistent with the 
ground state wavefunction. 

However, on other manifolds, such a transition may not change the number of 
loops. Consider, for example, a single closed loop (1,1) on a torus, winding once 
around the meridian and once around the longitude. This can transition to a loop 
with different relative windings (1,-1), as in the following hgure where opposite 
sides of the square are identihed 




77 


—)■ 

r 


Thus, the sign structure of the wavefunction on a torus obeys a more complicated 
rule. In fact, consistency of the transition rule on non-orientable surfaces can 
reduce the ground state degeneracy of the double semion model; see the next 
section. 


3. Some pathologies of DS and its potential generalizations 


For an arbitrary manifold S, let TC'(S) and DS'(S) denote the subspace of 
zero energy states. In this section we consider the degeneracy of this subspace for 
various manifolds. 

Closed non-orientable surfaces are diffeomorphic to: P, 2P, 3P, ..., where P is 
the projective plane, PP^, and the positive integer represents that number of Ps 
tubed (“connected sum”) together. 

Theorem 3.1. 

dim(DS(tP)) = 2*-^ 


Note. For orientable surfaces, DS and TC have equal dimension, but for closed 
non-orientable surfaces, ~ P 

We will prove theorem 3T later in this section. First we collect a few elementary 
topological facts and then prove some special cases as a warmup before proving 
the general result. 

Fact 3.2. If a closed surface is the boundary of a 3-manifold, 


dim her (Pi (S, Z 2 ) ^ Pi (M; Z 2 )) = - dim Pi (S; Z 2 ) 
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In particular, dimiJi (S; Z 2 ) is even. Equivalently the Euler characteristic x(S) 
is even. 

More generally, if is the boundary of a (2k+1)-dimensional manifold , 

> 1, then: 

dimker {Hk {N; Z^)) ^ H, (M; Z2) = ^ dimH, (iV; Z2) 

Again, this implies dimHk {N, Z 2 ) is even, and the Euler characteristic x{^) (base 
of Betti numbers with any coefficients) is also even. 

Since we only use Z 2 -coefficients, they are now dropped from the notation. 

Proof. Consider the long exact seqnences of the pair (M, N) in both cohomology 
and homology with the vertical isomorphism given by Poincare and Lefschetz 
dnalities: 


H^{M, d) 

L.D. 




L.D. 


hr^(s) - 

P.D. 




d) 

L.D. 


L.D. 

Hk{M,d) 


The universal coefficient theorem (with held coefficients F) provides a natu¬ 
ral isomorphism H*( ; F) Hom(iP*( ;F),F) implying that a and a above 
are Horn-dual. Thus Hom(coker a, Z 2 ) = kerd. In particular, dim(cokerQ;) = 
dim(kerQ;) = dim(ker 7 ), the later from the commutative square. By exact¬ 
ness of the homology sequence, coker a = image= keri*. Thus dim(keri*) = 
dim(ker 7 ). But using exactness once more, image i* = ker 7 so we conclude: 
dim(kerd) = dim(imaged). Since dimiPfc(E) = dim(keri*) -|-dim(imagei*), the 
fact is proved. □ 

Corollary 3.3. tP, t odd, do not bound any compact 3-manifold. For t even, tP 
does bound a compact 3-manifold. 

Proof. An elementary calculation shows Hi (tP) = Z\ and thus has odd dimension 
for t odd. 

Let K = 2P be the Klein bottle. iP is a S'^-bundle over with rehection 
monodromy. [S^ is the circle.) This monodromy clearly extends to a rehection of 
the disk D^, so K bounds a disk bundle over . This handles the case t = 2. In 
general, for t = 2j one may tube together j copies of this disk bundle to produce 
the desired 3-manifold. □ 


Fact 3.4. Any two embedded codimension = 1 submanifolds P and Q of a manifold 
which represent the same Z 2 -homology class a G H^-i {N; Z 2 ) are cobordant: 
There is a compact with = PWQ (and in fact W also maps into N 

extending the inclusions of P and Q). 
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We will apply this where P and Q are closed surfaces and N is a closed 3- 
manifold. 

Proof. Hd-i (iV, Z2) = {N, dN, Z2) by Lefschetz duality. is represented as 

homotopy classes of maps [N, RP°°] (with a relativization on dN to the base point 
if dN 7 ^ 0). We use dX to denote the boundary of any space X. The classifying 
space for ( ; Z2) is K {Z2, 1), the space with tti = Z2 and all higher homotopy 
groups vanishing; it is the inhnite projective space RP°^. Its characteristic class 
L G {RP°^-, Z 2 ) is dually represented by a codimension one RP°°~^ C RP°°. 
A class p e {N,dN; Z 2 ) is represented by a map fp : {N,dN) —)■ {RP°°,*); 
P. D.(p) = [f~^{RP)°°~^y The submanifolds P and Q arise as different but 
homotopic ways of making the representing map transverse to RP°°~^ C RP°°. 
W is produced by a relative application of transversality to the homotopy, F. That 
is, P = /-I Q = /-I and W = P-^ (PP~-i). □ 

Now consider DS(PP^). We hrst consider the question: Is there any ground 
state wave function (g.s.w.f.) ) G DS {RP^) with nonzero weight on the empty 

picture or blank multi-loop diagram. The answer is no, since starting with the 
empty picture, a loop may be swept over P with 3 local events, each of which 
produces a —1 phase difference between the multi-loops before and after each 
event. The product of these three signs is —1 showing '^(0) = —1/’(0) = 0. 



In words, we create a small circle in a Mobius band Xi, expand it around the 
band, and recouple it so that it becomes parallel to the boundary of Xi. Since 
P = Xi Ud disk, the recoupled loop now shrinks to a point across the disk and 
disappears. The sign is (—l)xT) as each “event” corresponds to a critical point of 
a Morse function on P. In the simplest case there are three. The crossing in the 
above hgure represents the half-twist in the Mobius band; the fourth image (out 
of six images in the sequence) represents the loop parallel to the boundary of the 
band. This boundary is drawn, after projecting the band into the plane, as a loop 
with a single self-crossing. 

Similar argument shows that for any DS g.s.w.f. '^tp(empty) = 0, f odd. The 
sweepout now has t + 2 events: the birth of a circle, t reconnections, and hnally 
the death of the circle. We have shown: 

Lemma 3.5. Fort odd, anyDS ground state wave function satisfies ^jJtp (empty) = 

0 . 


□ 

To gain some intuition, let us now consider a g.s.w.f. 'ipp( ) with ipp (RP^) = 1, 
which evaluates to one on the essential loop. Consider P as the unit disk D 
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modulo antipodal boundary identifications. There is an embedding of the Klein 
bottle K ^ P X corresponding to revolving RP^ (a diameter of D) by vr as 
one traverses the coordinate. 



Figure 3.1. 


Any 1-parameter history of multi-loops including births, deaths, and reconnec¬ 
tions beginning and ending with RP^ dehnes a closed surface S ^ P x S^. The 
second homology of P x obeys the Kiinneth formula: 

H2 (P X S^) = Pi(P) X Pi (5I) © H 2 {P) X Ho {S^) = Z2 © Z2 

Choosing the initial multi-loop to be RP^ or any other choice in that homology 
class hxes the hrst factor to be 1 G ^ 2 . The two remaining choices corresponding 
to H 2 {P) = Z 2 are exemplihed by the product history (nothing happens!) where 
S = X := T, a. torus, and the twisting history (Figure 3.1) where S = K. 

Any other choice of history S' is homologous in P x 5^ to either T or K. Since 
both T and K are cobordant to 0 (bound closed manifolds), and (by Fact 3.4) S' 
is cobordant to T or K, we see that S' is cobordant to 0. (Concretely, glue the 
cobordism from S' to T or iC to the null cobordism along their common boundary 
(T or K).) Thus by Fact 3.2, x{S') = even. Consequently the skein relations 
between multi-loops in the non-trivial homology class are consistent: an even 
number of (—l)s appear as we move from any essential multi-loop conhguration 
through others and finally back to itself. This means that there is a nonzero DS 
g.s.w.f., with 'ipp (PP^) = 1. With this warmup we can complete the proof of the 
theorem. 


Proof of \3.1\ Let S be a closed surface with hrst Z 2 -Betti number = b. Let the ele¬ 
ments Cl,..., 626 of Pi(S) span . There is a map 9, (Pi),..., "0 (P 2 O) 

which injects the Hilbert space DS(S) into , where Pj is any multi-loop repre¬ 
senting Cj. 

To determine dim(DS(S)), it is sufficient to characterize the e* for which {Ef) 
can be nonzero. 
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Claim 3.6. There is a DS g.s.w.f. 'ipi so that 'ipi (Ei) = 1 and 'ipi (Ej) = 0, i j, 
if and only if wi (cj) + x(S) = 0. wi is the first Stiefel-Whitney class. {0 {'ipi)} 
spans image 9. 

If S is orientable, note that wi (e*) + x{S) = 0 + 0 = 0 mod 2. 

Also note the claim implies the theorem since, for a non-orientable S = tP, 
Wi : Hi{'E) —)■ Z 2 is onto, so exactly | of the 2* first homology classes have 
(cj) = 0 and exactly I have Wi {ef) = 1. 


Proof of 3j_6. As in our warmup discussion of RP^, there are precisely two homol¬ 
ogy classes x in S x which realize a given Cj upon intersection with E x 1. 
One is represented by the constant history. Any surface S (history) in this class is 
cobordant to Ei x which bounds EiX D"^, so S is cobordant to 0, hence +(5') = 0 
mod 2. In this case, the skein relations merely say that 'ip {Ef) = (— {Ef) 
and pose no restriction. The other possibility is more interesting: E homologous 
to A’i X + E X 1. In this case we compute: 


(*) X{S) = w, {Ei) + x(E) mod 2. 

Ip {Ei) is forced to vanish exactly when the right-hand side is odd. Thus it suffices 
to establish 0 for any representative of the cobordism class of S. There is no 
loss of generality taking Ei to be a single circle, which we do. The natural choice 
is a “resolution” of the union {Ei x U (E x 1), a union of a torus and a copy 
of S along Ei X 1. This means that in each normal disk cross-section oi Ei x 1 
we resolve the crossing —|— of the two surfaces as either or _ yyg 
do this continuously all the way around x 1, and we will succeed precisely 
when Wi (cj) = 0, for in this case Ei G T, has a cylindrical neighborhood (as does 
Ei X 1 (Z Ei X S^). We will fail precisely when wi {Ei) = 1, i.e., Ei has a Mobius 
band neighborhood in E. In this case, near the contradictory point we see the 
two resolutions fitting together to form a “baseball” curve near the final point (see 
Figure |3.2l). 



Figure 3.2. 
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Thus, to compute the resolution, a disk 6 must be added bounding the baseball 
curve in the ball in Figure [3^ Using the additivity formula for Euler characteristic 

(**) + 

one readily checks 0 by showing that when wi (Ei) = 0, x(Res) = x(S) and 
when wi {Ei) = 1, the partial gluing along two arcs (not circles) contributes —2 
to X while 6 contributes +1, so in total x(Res) = x(S) ~ 1- ^ 

□ 

Note. As an example, on K any DS g.s.w.f. evaluates to zero on the two non- 
orientable 1-cycles, but is arbitrary on the oriented cycles. In general the restric¬ 
tion on a DS g.s.w.f. on tP is that it evaluate to zero on non-orientable cycles if 
and only if t is even and that it evaluate to zero on orientable cycles if t is odd. 

In subsequent sections we will identify classes of d-manifolds X for which GDSs 
on X have ground state degeneracy (gsd) = \Hi (X;Z 2 )|. In the remainder of 
this section we show a dramatic drop in ground state degeneracy if DS is gen¬ 
eralized to a theory of multi-loops of dimension one fluctuating within closed 
3-manifolds. This motivates our focus on fluctuating codimension one surfaces 
within d-dimensional manifolds. 

First consider a theory in which unframed multi-loops fluctuate in a 3-manifold 
with the zero mode ground state wave function (zgswf) experiencing a phase 
= — 1 with each Morse transition. Even locally, within a 3-ball, a closed cycle with 
one Morse transition exists: 



The local inconsistency: "0(circle) = —'(/’(circle) forces any zgswf to be identically 
zero on any M^. 

The explanation for this local inconsistency is that we neglected a normal fram¬ 
ing, which we now add. So we now assume all Morse transitions respect framings. 
The beginning and end pictures on line (3.1) now have different framings so there 
is no inconsistency. 

However, even in the context of framed multi-loops, we hnd that the Hilbert 
space V {M ) of this theory for any closed oriented 3-manifold M has dimension 
= 1, and below (Theorem 3.7) characterize the single nontrivial sector. 


Note. It is proved in [7] that for any modular tensor category (MTC), the skein 
space on a closed oriented 3-manifold has dimension one, so the dimension count 
in Theorem 3.7 is merely a special case, since “semions” form a MTC. 
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Given a closed 3-manifold, let E G Hi{M) = = horn -^ 2 ) 

denote the class which assigns to any closed surface embedded in M, S ^ M, its 
Euler class x(S) mod 2. 


Theorem 3.7. Let M he a closed oriented 3-manifold and let V{M) be the “Hilbert 
space” of framed multi-loops embedded in M with skein relation: “Morse = —1”. 
diml/(M) = 1, with E being the nontrivial sector. 


Proof. The argument closely parallels the proof of Theorem 3.1 
sketched. Let L be a loop or multi-loop in class E G Hi{M) 


so it is only 
Consider a “history” 


h beginning at L x 0 C M x 0 C M x [0,1] and ending at L x 1 C M x 1 C M x [0,1] 
and then identify boundaries to embed M x (0,1) M x [0, l]/m x 0 = m x 1 = 
M X S^. Let us compute the mod 2 Euler class of the closed history h (Z M x S^. 
[h] = [L] X [S*^] -I- a; X [pt.] G H 2 (M x S^), where x G H 2 {M). Using the data of 
the normal framing, one computes: 


x{h) =x{L X 5^) -I- x{x) + |L X n x| mod 2 
= 0-1- x(a;) + E ■ X mod 2 

The right-hand side is zero for all histories, i.e., for all x G H 2 {M), if and only 
if F = [L] = E, which is precisely the case that the sector survives in V{M). □ 

4. The theories we study: “fluctuating {d - 1) submanifold of 

A smooth closed manifold will be the home or “ambient space” of our gen¬ 
eralized double semion (GDS) theory. We consider a hxed combinatorial structure 
C, a generic cellulation, on X. Associated to C is a “crude Hilbert space” PL con¬ 
sisting of one qubit for each (d — l)-cell of C, spanned by (t= present, 1= absent). 
There is a Hamiltonian: 

H = H+ +Ha 

acting on Pi. We further write: 

H+= ^ H,,Hu= Y. 

d—2 cells e d-cells c 

H_|_ has a term for each (d —2)-cell P which is zero if an even number of (d — l)-cells 
meeting P are f (present) and 1 if an odd number are f. This is the same as the 
/d_|_ term for the generalized toric code. 

Ha has a term He for each d-cell c. This term is 

l-Oe 
2 ’ 

n Y 

1 ) cells f dc 


where 

(4.2) 


Or. = ± 


(d- 
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The operator Xf is the Paul i-X operator acting on the cell /, changing the state 
of the cell. The sign in (4.2) is —1 x —ix(tc) -^vhere x is Euler characteristic and 
tc is the codimension = 0 submanifold of dc consisting of the union of {d — 1) 
cells of dc which are labeled t in the state on which He acts. We will show in 
the appendix that this subset is indeed a smooth submanifold with corners. This 
is the reason that we chose a generic cellulation; it is analogous to the reason for 
dehning the double semion model in two dimensions on a hexagonal lattice, rather 
than a square lattice. 

The boundary of tc regarded as a submanifold of dc is the same as the boundary 
of j-c and is equal to tc H 4,c. For even d, the Euler characteristic of tc H j-c is 
equal to twice the Euler characteristic of tc, so — is equal to ix(tcn4.c)_ This 
makes the sign more reminiscent of the sign in the double semion model in two 
dimensions, where there is a factor of i for every leg leaving the hexagon. 


Explanation of generic cellulations. Generic cellulations may be dehned as those 
divisions of a smooth closed d-manifold X'^ into a union of smooth fc-cells, 0 < 
k < d, piecewise smooth on their boundaries, which obey the local combinatorics 
of the dual cells to a smooth triangulation of X. For example, on a surface 
(d = 2), a trivalent graph with contractible complementary regions determines a 
generic cellulation. In 3D the point singularities are cone over the 1-skeleton of 
a tetrahedron, etc. A second construction which works in any dimension is to 
start with not a smooth triangulation but rather Riemannian metric. Then any 
sufficiently dense set S of points, if perturbed to be generic, will determine Voronoi 
cells (the d-cell interiors consist of all points closest to some s G S') which dehne a 
generic cellulation. We work with these structures rather than some hxed lattice 
for two reasons: 

(1) every union U of (d — 1) cells whose Z 2 -boundary vanishes every (d —2) 
cell meets an even number of (d — 1) cells of U) is a submanifold. 

(2) heritability: The intersection pattern induced on the boundary of any cell 
of a generic cellulation is also a generic cellulation in a lower dimension. 
For more details, see appendix. 

Lemma 4.1. All terms of pairwise commute. All terms in id+ commute 
with all terms in Hq. All terms in idn pairwise commute when restricted to the 
eigenspace of with vanishing eigenvalue. 

Proof. The terms in pairwise commute since they are all diagonal in the prod¬ 
uct basis I 'O, I J,). We assume the reader is familiar with the analogous verification 
that the terms of Kitaev’s |2] toric code (TC) commute. As with TC, adjacent 
-|- and □ operators share two anticommuting factors and therefore commute. The 
local sign is irrelevant. 

Consider a pair of terms in . They trivially commute unless dci and 

dc 2 share some (d — 1) cell /. Consider a state T in the product basis on which 
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vanishes. We show that = 0c20ci^. The states and Oc20ciT 

are both product states and have the same spin conhguration. However, we must 
check that they have the same sign. Suppose hrst that / is labeled in T, so that 
Oc 2 changes the set of t spins on ci to tci U/. So, the sign for Oc^Oc 2 is equal to 
_]^x(tciU/) _ ]^x(tc 2 )_ similar computation for Oc 2 i^ci gives — 1 ^(H 2 W) _ ]^x(tci)_ 
So we must check that 

(4.3) - = -1v(H 2U/) _ ix(tci)_ 

By additivity formula of Euler characteristic, y (tci U/) = y (tci) + x(/) — x(tci 
n/). So, we must check that 

(4.4) x(tci n/) + x(tc 2 n/) = 0 mod 2. 

However, by assumption that vanishes on T, '1'^ H/ =tc 2 hi/. If instead / 
is labeled / in we must check that —— l^(l''= 2 ) = _xv(tc 2 \/) _ xx(t':i). 
A similar calculation using additivity formula for Euler characteristic shows this 
case also. □ 


Lemma 4.2. All terms of H are projectors. 


Proof. This is obvious for terms of i/+. For He one needs to check that the operator 
Oc is order two. To do this, one needs to check that the signs agree: 

(4.5) (—1)^'" = (—1)^'*"^ mod 2 


where /c is the complement in dc of /c- But this follows immediately from three 
facts: 

(1) dc is a (d — l)-sphere xi^c) = 

(2) the additivity formula: 


0 d even 
2 d odd ’ 


X (tc) + X (ic) - X (tc n le) = xi.de) = 0 mod 2, 


and 

(3) Since (/c h ic) bounds (it in fact bounds both fc and /c), X (tc h /c) = 0 
mod 2 since the Euler characteristic of any bounding manifold is even. 


The upshot is (4.3). 


□ 


As in the double semion case in two dimensions, it is possible to dehne a Hamil¬ 
tonian which is a sum of commuting projectors by projecting the operators He 
into the zero eigenspace of the terms He for e a face of c, setting Hu = 
with 

(4.6) 

ffr”=( n i-i4.)ffo( n ^-H,). 

d — 2 cells e, s.t. e is a face of c d — 2 cells e, s.t. e is a face of c 


Then, 
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Lemma 4.3. All terms pairwise commute and commute with all terms He- 


Proof. The proof is the same as lemma 4.1 
vanishing of He for He a face of Ci, C 2 . 


noting that that proof only nsed the 

□ 


Explanation of the sign of He- We shonld think of He as indncing flnctnations 
between codimension one snbmanifolds. The most basic possible flnctnations are 
Morse transitions which near the transition point, by dehnition assnme the follow¬ 
ing local form, for index 0 < k < d — 1: 

2 2 2 I 2 I I 2 , 

—Xi — X2 — ■ ■ ■ — Xf, + x,._^_i -t- • • • Xei_i — —e ^ 

-xl-xl - xl + 4+1 • • • -f 4-1 = e, e > 0. 

For example, for the nsnal donble semions these are the transitions: 

0^O(fc = o) 

O ^ 0 (i = 2) 

Form the usual Skein relations, each of which is associated with a phase = — 1. 
We wonld like to generalize this to higher dimensions by declaring a phase = — 1 
for each Morse transition. 

Unfortnnately, the transition t :=tc~^4'c is not in general Morse bnt in general 
mnch more complicated. However, a simple argnment (below) enables ns to write 
t as a composition of (y (tc) + l)mod 2 Moi’se transitions. This decomposition of t 
is not nniqne bnt the parity of the nnmber of Morse factors is. 

Consider a Morse fnnction: / : c —>■ [0,1] with /“^(O) = (tc)~, = (ic)”, 

the indicating deletion of a thin collar. Let / be free of critical points on 

Sc\((td-U(«-). 

As we cross a critical level s of /, y (/ ^ (0,s+]) changes by ±1 according to 
index(s) = odd or even, respectively. Since y(c) = y (/“^[0,1]) = 1, the nnmber 
(mod 2) of Morse transitions of / is y (tc) + 1 = X (/~^(0)) + 1- 

Observation 4.4 (“Morse = —1”). The Hamiltonian H enforces on any zgswffj, 
the skein rule Morse = —1. That is, if E and E' are cellular d — 1 submanifolds 
of M which differ by k Morse transitions, '4){E) = {—lYf {E'). 


5. Some properties of CDS 

Let X'^ be the ambient manifold of dimension d. As before, onr model is reg- 
nlarized by choosing a generic cellnlation C of X. In this section and Section 
cntting and glning constrnctions are nsed whose motivation is from the contin- 
nnm. It may be necessary to snbdivide C to carry ont these constrnctions. Thus 
all statements should be interpreted as holding for some rehnement C' of C. 
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The first property we observe regards the form of zero energy ground state wave 
functions (zgswfs). Because of the term all zgswfs are superpositions of (cel¬ 
lular) submanifolds (of dimension d — 1, and not necessarily oriented). Because 
Hu does not fluctuate between different classes of Hd-\ (X; Z 2 ), all zgswfs are su¬ 
perposition of zgswfs each with support confined to a fixed class x € Hd-i (X, Z^. 
The next lemma implies that there is at most one zgswf up to phase for any class 
X G Hd-i (X, Z 2 ). However, as we saw when d = 2 on non-orientable surfaces, 
it is not necessarily the case that each possible sector x G Hd-i (X;Z 2 ) actually 
admits a zgswf. 

Lemma 5.1. Let a and {3 be two C-cellular (d — l)-Z 2 -cycles belonging to the 
same class x G Hd-i (X; Z 2 ). There is a finite set {ci,..., Cj} of d-cells of C so 
that a — fi = d {ci U ■■■ U Cj). 

Proof. This is directly from the definition of cellular homology. No multiplicities 
or orientation need be considered since the coefficients are Z 2 . □ 

Corollary 5.2. Let be the Hilbert space spanned by the elements ei,... ,626 
of Hd-i{X), and choose representative cellular cycles Ei for Ci, 1 < i < 2^. The 
map 9, (Ei ),..., "0 {E 2 b)), is an injection of the zero-modes GDS(X) —)■ 

C20 □ 

In Section]^ we discussed “missing sectors” when X is a non-orientable surface. 
In higher dimensions, there may be “missing sectors,” 9 has cokernel, even when 
X is orientable. An extreme example is X = CP^. 


Example. Let the ambient manifold X be CP^, complex projective 2-space, d = 
dim (CP^) = 4. P 3 (CP^, Z 2 ) = 0, so there is only one possible sector of zgswfs in 
CDS, the trivial sector. But, in fact, GDS (CP^) = 0. Just like PP^, CP^ has a 
Morse function with three critical points 


(^ 1 , Z 2 , Zfi) —)■ 


\Z2? + 2\Z. 


X1I2 + IZ2P + IZ3 


Thus, as in Sectionj^ 0(0) = —0(0) = 0, so GDS (CP^) = 0. Of course H will have 
some least energy state but we do not know that it has topological significance. 
We also do not know its entanglement properties. If we work with a Hamiltonian 
that is a sum of commuting projectors using terms then the ground state 

is an eigenstate of every term separately. Hence, this ground state will be a zero 
energy eigenstate of the Hamiltonian on some “punctured cellulation”, obtained 
by removing the terms Hc,He which have expectation value -|-1. We now give 
precise conditions under which 9 is an isomorphism: 


Theorem 5.3. If d is odd, 9 is an isomorphism, GDS(X) = Moreover, in 
each sector x G Hd-i{X), an explicit zgswf is given as 0x(P) = 
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For d even and x E Hd-i(X), there is an zgswf 'ipx{E) 7 ^ 0, [E] = x E Hd-\{X), 
if and only if x{X) + e(x) = 0 mod 2, where e{x) is defined below. 

In general, there is no local formul^ for the zgswfs; however, if X is a Z 2 - 
homology d-sphere^ the formula for the sign of ijjx on E is ± s{E) 

also defined below. 


Proof. Consider d odd. In this parity, the codimension one cellnlar snbmanifold 
E is even dimensional. Morse transitions of even dimensional manifolds always 
change y by ±2. The transition removes an x and replaces it with 

Di+i X 


x{sn 


0 p odd 
2 p even ’ 


X (S> X 0“-'-^) = X (s^) =2 X (8“-’-^] = X (8“-’-^ X W+']. 
The claim now follows from the glning formula for Euler characteristic: 


(5.1) x{AUB)=x{A) + x{B)-x{AnB) 

Thus every Morse transition M ^ N changes iA^) by a factor of —1 (which is 
consistent with the fluctuation of Hq). 

Now consider d = even. Let E be a cellular hypersurface of X representing x. 
Let u be the normal bundle of E in X and L C E be ad — 2 dimensional embedded 
submanifold Poincare dual to the hrst Stiefel-Whitney class Wi{u) E H^{E), of the 
tangent bundle to E. 


Definition. e{x) = x{E) mod 2. 
Lemma 5.4. e{x) is well-defined. 


Proof. According to (3.4), any other hypersurface E' representing x is bordant to 


E via V C X X I; d+V = E', d-V = E. Furthermore, denoting the normal bundle 
to E' (in X X 1) by u', for any E P.D. to wi (z/') in E' there is a bordism (it is 
the P.D. to Wi {vv^xxi)) L <Z V from L to L'; dL = L]JL'. Thus L while not 
uniquely dehned is well-dehned up to bordism. 

From Fact |3.2 , if a closed manifold L of dimension 2k is a boundary, then bk{L) 
is even (and dim(kerfc) = An immediate corollary is that if is a closed 

manifold which bounds, then x{L) = even, and in particular, if is bordant to 


then x(L) = y (E) mod 2. 


□ 


^Even when d = 2, it is only true on the 2-sphere that the signs zgswfs are (_i)tt(components)^ 
Recall the example on the torus E, where the curves (1,1) and (1, —1) are in the same Z 2 -sector, 
but 1)) = —'b((l, —1)) for any zgswf tp. 

At is sufficient to assume Hd-i{X) = 0 = Hd{X). 








DOUBLE SEMIONS IN ARBITRARY DIMENSION 


16 


Now 'ipx{E) can be normalized to 1 unl ess there is some odd length seqnence of 
Morse transitions in X starting with E and retnrning to E. That is, nnless there 
is a snbmanifold W <Z X x with hh fl (X x 0) = and x{^) = odd. 

Bnt the Kiinneth theorem tells us there are exactly two possible homology classes 
(and therefore bordism classes) for W. They are x x [S*^] G Hd-i{X) ® Hi (S*^) C 
Hd{X X X X [^1] + [X] X * e Hd-i{X) ® (^1) © Hd{X) ® Hi {S^) = 

Hd (X X S^). Again, to determine the parity of x(hT), it is sufficient to study any 
representative of each homology class, since all representatives are bordant. In the 
hrst case, E x is & representative with y x = 0, so we may dismiss the 
hrst case entirely. 

The second case is more interesting and fully analogous to the discussion at 
the close of Section The class x x [S*^] + [X] x * is represented by the union 
£' X 5^ U X X 0. This union is not a submanifold but instead is singular along 
£'x0 = i?xS'^nXx0. If the normal structure, a “ —|— ” bundle over E x 0, is 
trivial (i.e., a product bundle), the singularity can be resolved (replace —|— with A^) 
with y(resolution) = x{E x S^) + y(X x 0) = y(X x 0). While the S'^-direction 
(“vertical”) segment of “ —|— ” is trivial, the horizontal segment is reflected along L; 
thus the resolution is as previously pictured in Figure 3.2 but now in a parameter 
family over L. Where previously a single resolution as in Figure [3^ accounted for 
an additional 1 mod 2, now, in its parametric form, it contributes y(T). Thus 
x(W) = y(X) + e(a;). The sector x disappears from GDS(X) precisely when y(fF) 
is odd. 

Thus completes the proof of 5.3 except for the statement about Z 2 -homology 
spheres which is discussed in Remark |5.6[ □ 


To make further progress, we need the notion of “Kervaire semi-characteristic.” 
We assume d = 2k + 2 is even. Kervaire dehned the semi-characteristic of a 
closed odd dimensional manifold 

k 

(5.2) s{M) = J2HM), 

6=0 

where bi is the Betti number of M. Kervaire considered real Betti numbers but 
it is more consistent with our approach to take Z 2 -Betti numbers: 

(5.3) b, = dim {H,{M-,Z 2 )) 

and consider s{M) G Z 2 , the mod 2 reduction only. 

The next lemma says that s fluctuates in a manner similar to the sign 6c of He. 


Lemma 5.5. Suppose undergoes a Morse transition to N of index < k. 

Then s{N) = s{M) + 1. If index = k, this formula still holds if either 

(1) the vanishing (k — 1 )-cycle 6 at the transition is nontrivial in Hk_i (M; Z 2 ) 
or 
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(2) [5] = 0 but the bordism W associated to the transition has vanishing inter¬ 
section pairing Hk+i (W] Z 2 ) x i/fc+i [W] Z 2 ) —)■ ^ 2 - 

The intersection pairing clearly vanishes on Hk+i {W ; Z 2 ) ximage {Hk+i (M; Z 2 )) 
so it is sufficient to consider the intersection pairing on the quotient 

W := Hk+i (IT; Z 2 ) / image (M; Z 2 ) = Z 2 


(5.4) 


W X W ^ Zj. 


Since we are working with Z 2 coefficients, the vanishing of the pairing (5.4) is a 
nontrivial assumption for W of all even dimensions. 


Proof o/ |5.5| From here forward, Z 2 coefficients are implicit, although most asser¬ 
tions and diagrams have analogs over other helds and sometimes rings. 


j + 1-handle 


X T)2fc-i+i^ 
dDi+^ X T)2fc-i+i^ 


^2k+l 





Figure 5.1. Curved dotted line is called “co-core” and hor¬ 
izontal dotted line /} 2 fe-i- 2 -o+i) jg called “core”. Solid dots indicate 


5.1 


The Morse transition of index j -|- 1 creates a cobordism W pictured in Figure 
The associate “braided” exact sequences of pairs is displayed below. 


0 0 

II? II? 


Hy+2{W,M)^Hy+,{M) 

if,+1 (IF, TV) 


Hj(W,N) 

X 

X 

X 

X 


(IF) 

H,{W) 

a ^ 

Hy^2{W,N)^H,^,{N) 

X 

ii,+i(lF,M) 


N, 

Hj(W,N) 


83 


II? 

^2 


H,{W,M) = 



*=j + l 
* ^ J + 1 


and H,{W,N) = 



* = 2k — j 1 

* 2k — j -\-l 
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First assume j < k making H^,(W,N) = 0, * < j + 2. Then Hj{N) = Hj(W). 
Suppose dj+i = 0, then also Hj{M) = Hj{W), so Hj{N) = Hj{M) and Hj+i{N) = 
Hj^i{W) = Hj+i{M) © Z 2 . On the other hand, if dj+i ^ 0 then Hj{M) = 
Hj{N) © Z 2 and Hj+i{M) = Clearly for i < fc and i < j or j + 1, 

Hi{M) = HiiW) = Hi{N). Thus in either case, dj+i = 0 or dj^i 7 ^ 0, the 
semi-characteristic s E Z 2 satisfies s{N) = s{M) + 1. 

Now consider the more delicate case j = k. 

Hk{N) ^0 

X 

Hk{W) 

X 

Hk{M) 0 

Again consider cases. For d, d' zero or nonzero, we have the following table for 
s{N)-s{M): 


a' = 0 

ax 0 

The * indicates that this case does not actually occur. It is forbidden by Poincare 
duality. If 5(1) = 1/ 7^ 0 G iffc(M), y will be Poincare dual to some y G Hk+i{M); 
y ■ y = 1. But a punctured copy of ^ is a null homology for a'(l) as sketched in 
Figure |0 


5 = 0 5^0 


0 

1 

1 

* 


0 - ^ Z 2 

X X" 

Hk+i{W) 

^ X 

0- ^ Hk^i{N) Z2 


handle 



2/punctured 


Figure 5.2. 


The (0, 0) case occurs precisely in the forbidden case where the cobordism W has 
a (fc+l)-cycle with odd self-intersection, or from line ( |5.4 ) Wx W —)■ Z2 is nonzero. 
One way to compute this is to take the obvious dual disks (called “core” and “co¬ 
core”) X 0 and 0 x in the A; + 1-handle x x of 

W relative M and close them off by chains, one in M exhibiting a null homology 
for 5(1) and one in N exhibiting a null homology for 5'(1). This is illustrated in 
Figure |5^ where W is drawn as a punctured Mobius band. 

□ 
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W 



single transverse self intersection point 
of the nontrivial element q G W. 

Both dark loops represent the class q. 


Figure 5.3. 


Remark 5.6. This remark completes the proof of Theorem 5^: the statement 
about manifolds X with H^{X) = (*5''^)- When Hd-i{X) = 0, each elemen¬ 

tary bordism W across a plaquette (i.e., d-cell) will necessarily be orientable and 
hence (unlike the bordism in between (1,1) and (1 ,—l)j be represented as a 
codimension zero submanifold of X , W <Z X. This and Hd{X) = 0 implies that 
even when Morse transitions are of index k, d = 2k -\- 2, condition of (5.5) 
will hold, whether or not [5] = 0. So we conclude that all Morse transitions of 
E C X change the parity of s{E). Furthermore, in the previous part of the proof 
we established that x{W) = even in the present case since x{^) = X {S‘^) = 2 and 
e(0) = 0. These two facts imply that (—1)**^®^ is a consistent assignment of signs 
for a nontrivial (and unique up to scale) zgswf in GDS(X). This completes the last 
portion of Theorem 5.3 Note that the proof only requires Hd-i{X) = 0 = Hd{X), 
a weaker assumption than assuming X is a Z 2 -homology sphere. 


Remark 5.7. Ford odd. Theorem |5.3| gives an explicit form of the ground state 
wavefunction in each sector as fJx{E) = We claim that the unitary U that 

multiplies each basis vector by a phase can be written as a local quantum 
circuit, i.e., with a bounded depth of the circuit and with the gates of the circuit 
supported on sets of bounded diameter, with these bounds depending only on local 
geometry. Thus, there is a local quantum circuit that maps the ground state sub¬ 
space of the GDS for d odd onto the ground state subspace of a GTG. Further, 
since every Morse transition for odd d changes x by ±2, if we restrict to the ground 
state subspace of then the unitary U conjugates the operator Hu of the GDS 
to the operator Hu of the GTG (we leave it as an open problem whether there is 
a relationship between the theories if we consider the space of states outside the 
ground state subspace of H+). 

To construct the local quantum circuit, write 

d-l 

(5.5) n ui.c), 

j=0 j-cells c 

where the second product is over all j-cells c, and U{j,c) is a unitary that that 
multiplies basis vectors by phase equal to either 1 or with the phase being 1 
if the given cell c is not in E and ±i if c is in E (the phase being +i if j is 
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even and —i if j is odd). Each U{j,c) is supported only on the {d — 1) cells that 
intersect the given j-cell c. Those unitaries are the gates in the guantum circuit. 
Unitaries supported on disjoint sets can be performed in the same round, and in a 
bounded number of rounds all unitaries can be performed, assuming the cellulation 
has bounded local geometry, i.e., assuming there is a bound on the number of {j + 1) 
cells meeting any j cell and also a bound on the number of j-cells in the boundary 
of any (j + 1) cell and a bound on the diameter of any cell. 

6. Generalized loop, or “balloon” operators 

Let us recall the Wilson loop “W” operators |S] for the usual TC and DS 
theories. In the TC case, one merely specihes a loop /, the support of hfo, and 
then Wi acts as NOT on the bonds of 1. In the case of DS, this dehnition of Wi is 
augmented by a sign rule (e.g., line 44 of [5]) which is necessary to preserve ground 
states. The rule stated in [5] for a loop I or the 2-sphere S'^ may be restated in 
topological language as the sign rule below. 

Let a be the multi-loop of f bonds. The sign for the action of Wi on 'ijj{a) is: 


Again s is Z 2 -semi-characteristic, which is 1 for a simple closed curve 1. {I fl a) 
consists of line segments (or all of / if / C a), so d{l fl a) is an even number of 
points and x nierely counts them. These points are naturally divided into pairs 
(“0-spheres”) as the boundaries of the component arcs of a \ Z. For a collection of 
0 -spheres in a circle I, there is a total mod 2 linking number which we denote by 
“link,” which can be computed, for example, by making I bound a disk A, then 
spanning all the 0-spheres by simple arcs in A, and hnally counting (mod 2) the 
number of pairwise intersections of the arcs. An equivalent dehnition is to rehect 
all component arcs in a \ / which lie outside A into the interior of A, and then 


count pairwise intersections. An example is given in Figure 6.1 


I (round) /a 



a (in bold) 


Wi 



a 


sign = (-1) (-1^) = 1 


Figure 6.1. 

Our task in this section is to dehne the correct signs for the analogous Wilson 
“balloon” operator. To obtain explicit formulas, we introduce restrictions on the 
topology of X as necessary. These operators are equal to ± NOT on an imbedded 
codimension = 1 C-cellular submanifold C A'^, up to a sign. To be consistent. 
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the sign ± must be set so as to carry zgswfs to zgswfs. Keeping similar notation, 
a is the {d — 1)-dimensional cellular submanifold supporting a basis vector and L 
the Wilson balloon. 

Note that in addition to these Wilson loop operators Wi, the operator given by 
the product of Pauli Z around any closed loop on the dual lattice also commutes 
with the Hamiltonian of both the TC and DS models. This operator (which we 
will call a dual Wilson loop) straightfowardly generalizes to the GTC and GDS. 
For the GDS, given any closed loop on the dual lattice, the product of Pauli Z 
over all {d — l)-cells that intersect that loop is an operator that commutes with 
the Hamiltonian. If the closed loop is nullhomologous, then this operator acts as 
the identity on the ground state subspace. 


Theorem 6.1. Let d = 2k + 2, k > 0. As in 5.3, we require that X'^ is a Z^- 
homology sphere]^ Without this assumption, a local description of the balloon op¬ 
erator signs may not be possible. The consistent sign for Wilson balloon operators 


is: 


^^x(9(Lno))^ 


For the case of odd dimensional d, let d = 2k + 1. 
restriction on topology), the consistent sign is: 

jx{L) (^_l^xiLna) 


Then for all closed X'^ (no 


Note. The formula in odd dimensions is quite elementary, involving only Euler 
characteristic. In even dimensions, the formula is curiously simpler than the clas¬ 
sical d = 2 case of DS. If one traces the reason back through the proof, the 
anomaly for d = 2 is that d{L fl a) is then zero dimensional, i.e. points. In this 
very low dimension, the middle dimensional (zero) homology of S'® x S'® a product 
of spheres has rank 4 not 2 and the usual analysis of “hyperbolic pairs,” forms 


0 1 
±1 0 


or over Z 2 simply 


0 1 
1 0 


, does not apply. 


Proof of \6.1\ Consider hrst d odd. Write L = AU B, and a = B U C, M = 
dB = dC, as shown in Figure and apply the additivity formula (5.1) for 
characteristic. 


dA = 
Euler 



Figure 6.2. 

®It is sufficient to assume = 0 = Hd{X) and that the Stiefel-Whitney classes Wi of 

the tangent bundle t{X) vanish for z < |. 
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( 6 . 2 ) 

X(A U C) - x{B UC)= x(A) + - X{M) - x{B) - x{C) + x{M) = 

x(^) - x{.B) = xiA) + x{B) - 2xiB) = x{A U B) + xiM) - 2x{B) = 

x{L) + x{,M) -2x{Lf\a). 


Since M is a closed odd dimensional manifold, x{^) = 0- Thus the change in sign 
of a conhguration in the gswf upon applying Wl is as claimed. 

Now consider d = 2k + 2, even. To handle this case, we need a rather powerful 
algebraic lemma 6^ which, since we are using Z 2 -coefiicients, may be stated uni¬ 
formly regardless of whether k is even or odd. Lemma 6^ will be applied to the 
three kernels 


ker {Hk{M) ^ Hk{A )), ker ^ Hk{B )), and ker ^ Hk{C)) 

which, in a slight abuse of notation, are also denoted A, B, and C, respectively. 

Lemma 6.2. If the Stiefel-Whitney classes of the tangent bundle Wi (r 
vanish for i < /c = (| — l), then the intersection form on (M^^; Z) is isomor¬ 
phic to a direct sum Q = 0^^^ 

Proof. The total Stiefel-Whitney class of t{M), w = l-|-r(;i-|-r (;2 + -- -+ W 2 k, 
satisfies w = Sq(n) where v = l + vi+V 2 -I--- + V 2 k is the total Wu class and Sq 
is the total Steenrod square Sq = 1 -|- Sq^ -|- Sq^ Inductively we hnd: 



Wl = Vi 


W2 = V2 -I ViU Vi 

W 3 = Vs -\- Bockstein(n 2 ) 

W/i = V 4 + Bockstein(n 3 ) -f- ^2 U ^2 


so Wi = 0 , i < k, implies Vi = 0 , i < k. 

We have assumed that Wi{T{X)) = 0, i < k. However, t{X)\m = r(M) © 
and vm^x is trivialized by the cross-sections (inward normal to A, inward normal 
to B). Thus, for all i, WiT{M) = inc* Wj(r(X)), implying that Wi{T{M)) = 0, 
i < k. Hence the Wu class of M, = 0 G vanishes. 

Sq’^{x) = Vk U X for X e H^{M) so: 

xU X = Vk^I X, 

so the vanishing of Vk implies that for all x G x U x = 0. It is now an 

elementary exercise in duality to put the intersection form of M in the claimed 
coordinates. □ 
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In Lemma | 6 . 3 
a field. 


Z 2 is called F 2 , the field of two elements, to emphasize that it is 


Lemma 6.3. Consider a space , with coordinates labelled 1,2, Let Q be 

the form 


(6.3) 


Q 


/o 1 
1 0 

0 

1 

V 


\ 

1 

0 

•••/ 


Let A,B,C be three maximal subspaces (i.e., j-dimensional subspaces) such that 
Q vanishes on each subspace. Assume also that for any triple of vectors a E A,h E 
B,c E C such that a + b + c = 0 we have 


(6.4) (a,Q 6 ) = 0. 

Then 


(6.5) dim(y4 n S) + dim(i? 0(7) + dim(C' O A) = j mod 2. 


Since we will give an inductive proof and we will want to show that this last 
condition Eq. (6.4) holds under the induction, we will be refering to it below; 
so, we give it a name, referring to it as the “no twist” condition. Before giving 
the proof, we remark that the no twist condition is symmetric under interchange 
of the subspaces A, B, C, because ( 6 , Qc) = ( 6 , Qa) + (b, Qh) = 0 and similarly 
(a, Qc) = 0. 


Proof. Assume that dim(A n B) 7 ^ 0. In this case we will give an inductive 
proof, reducing the problem to a problem with j reduced by 1. Apply a basis 
transformation so that A and B both contain the vector vi = (1, 0, 0,0,...). Then, 
every vector in A, B vanishes on the second entry. So, A can be written as the 
span of Vi and A' and B can be written as the span of Vi and B', where A' and B' 
are j — 1-dimensional subspaces with Q vanishing on A! and B' and every vector 
in A' or B' vanishing on the hrst two entries. Suppose that vi is also in C or that 
V 2 = (0,1,0,0,...) is in C. Then, C can also be written as the span of vi and 
C or V 2 and C", with C vanishing on the hrst two entries. Then, the subspaces 
A', 5', C are maximal subspaces of a space with Q vanishing on A', S', C 

and the no twist condition holding for A',B', C because A', B', C are subspaces 
of A, B,C respectively. Further, 

, ^ dim(A n S) -I- dim(S fl C) -|- dim(C fl A) = 

^ ' dim(A' n B') + dim(S' fl C') + dim(C" fl A') -|- 1 mod 2 , 

as desired. 

So, suppose that Vi and V 2 are not in C. Then, it cannot be the case that every 
vector in C vanishes on the hrst coordinate (else V 2 must be in C) or that every 
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vector in C vanishes on the second coordinate (else vi must be in C). So, one 
of two cases hold. Either C contains a vector w = (1, 0, tcs, tC 4 , ...) and a vector 
w' = ( 0 , IjtCg, w; 4 ,...) or C contains a vector x = (1, l,a; 3 ,a; 4 ,...) but does not 
contain any vector of the form w or w'. 

Suppose the hrst case holds, that C contains such a w and w'. So, C can be 
written as the span of w,w' and some subspace Cq that vanishes on the hrst and 
second coordinates. Let y be the vector (0, 0, ^ 3 ,^ 4 ,...). Dehne C to be the span 
of Cq and y. Note that y is not in Cq by the assumption that Vi is not in C so C 
has dimension j — 1. Note that dim( 44 nC) = dim(y4nspan(r(;, Cq)) = dim( 74 'nC") 
and similarly dim(i?' fl C) = dim(i? fl C). So, this again gives three subspaces 


A',B',C' for which Eq. ( 6 . 6 ) holds and the no twist condition holds so again we 
reduce j by 1 . 

Suppose instead the second case holds, that C contains a vector 

X = (1,1,0:3,0:4,...) 

but does not contain any vector of the form w or w'. Dehne C to be the projection 
of C onto the space on which the hrst two coordinates vanish; i.e., C is the space of 
vectors (0, 0, 0 : 3 , ^ 4 ,...) such that (0, 0, T 3 , 4 : 4 ,...) is in C or ( 1 , 1 , 4 : 3 , 4 : 4 ,...) is in C. 
Being a projection, C has dimension j or j — 1. The form Q vanishes on C so C 
has dimension j — 1 . So, if a vector ( 1 , 1 , 0 : 3 , 0 : 4 ,...) is in C then (0, 0, 0 : 3 , 0 : 4 ,...) 
is in C. So, (1,1, 0, 0,...) is in C so C is the span of (1,1, 0, 0,...) and C", so 


dim( 74 ' n C) = dim(y4 fl C) and dim(i?' fl C) = dim(i? fl C) so Eq. ( | 6 . 6 [ ) holds. 
Further, the no twist condition holds for A', B',C'. 

So, we succeed in reducing j by 1 whenever dim( 74 ni?) 7 ^ 0 . A similar reduction 
can be performed whenever dim(A fl C) 7 ^ 0 or dim(i? fl C) 7 ^ 0 so we can assume 
that Ani? = i?nC = CnA = 0. We now show that this implies that j is even. 

Make a basis transformation leaving Q invariant so that A is the space of vectors 
vanishing on the even coordinates. Introduce notation, writing a vector v G as 
V = (t; y) where x, y are j-dimensional vectors to indicate that the odd coordinates 
of V are the entries of x while the even coordinates are the entries of y. Parametrize 
B as the space of vectors [Ku] Lu) for two matrices K, L where 44 is a j-dimensional 
vector. Since B is j-dimensional and A A B = 0, L is non-singular. So, we 
can instead parametrize B as the space of vectors {KL~^y,y) or (My^y) where 
M = KL~^. Similarly, we can parametrize C as the space of vectors {Ny]y). 
The matrices M, N are symmetric because Q vanishes on B, C. We now use the 
no twist condition. We have ((M -|- N)y,0) + {My,y) + {Ny,y) = 0, and so 
{My, y) -I- {Ny, 4/) = 0 for all y or equivalently {y, {M + N)y) = 0. So, M -|- 
vanishes on the diagonal. However, M+N must be nonsingular (otherwise, BAC 7 ^ 
0). So, i must be even. 

The fact that j must be even follows from the fact that any symmetric, non¬ 
singular matrix over F 2 which vanishes on the diagonal must be even dimensional. 
To show this, consider any such matrix F. If G is a non-singular matrix, then 
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G^FG is still a symmetric non-singular matrix which vanishes on the diagonal. 
We use two different choices of G, either a permutation matrix or a matrix which 
is 1 on the diagonal and has one nonzero entry of the diagonal, to simplify the 
matrix F through a sequence of steps similar to Gaussian elimination; the second 
choice of G gives a “row and column” operation. F must have a nonzero entry on 
the hrst row (otherwise it is singular), so use permutations to make the upper right 
corner of F nonzero. Then, use “row and column” operations to make all entries in 
the last column vanish, except for the entry in the upper right corner. Repeat this 
procedure on the second row, and so on, continuing until all the nonzero entries 
are on the diagonal from top right to bottom left corner. This diagonal intersects 
the main diagonal if j is odd (and so the matrix must be singular), so j must be 
even. □ 

Again, the notation of the lemma is designed to evoke Figure [6^ with a minimal 
abuse of notation. We wrote A, B, and G for the kernels of the three maps: 

Hu{M) ^ Hk{A) 


Hk{M) ^ Hk{B) 

Hk{M)^Hk{G). 

As we now show, these kernels indeed have the “no twist” (nt) property: 


(nt) 

Lemma 6.4. 


a + b + c = 0 {a G A, bGB,cGG) implies (a, Qb) = 0. 


Under the assumption Wi(r(X)) = 0,i<k + l = ^, (nt) holds 


Proof. First recall from the proof of Lemma 6.2 that Wi{T{X)) = 0, 
implies that two Wu classes vanish, Vk+i C and Vk G Fl^{M). The hrst 

vanishing is equivalent to x ■ x = 0 for all x G iffc+i(X) and the second vanishing 
is equivalent to y ■ ?/ = 0 for all ?/ G Hk{M). 

Assuming {a,Qb) = 1, we may build a cycle x G Hj{X) with x ■ x = 1, a 
contradiction. The construction is to form V representing x by gluing together 
three coboundaries: V = U 6 ^ U where the superscript is the ambient space 
of the coboundary, and d = p). 

The plan is to compute x • x taking V and a copy V perturbed in X to be in 
general position with respect to V and counting intersection points. The answer 


IS X • X = 




= (a, Qb), but from our assumptions on X, x • x = 0, implying 

(a, Qb) = 0. 

Although the submanifolds A, B, and G play symmetric roles, it is helpful for 
the computation to break symmetry and consider aUb(ZL-.= AVJB. We may 
obtain a second copy a' U 6' of a U 6 by perturbing this fc-cycle toward A using the 
normal direction of M C A. Let , 6 ^', and be formed by extending this 
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this crossing leads to an 
intersection point of a; ■ a; 


Figure 6.3. 

perturbation over the three Ss. There is a Z 2 -linking number dehned: 

/ (a U 6, a' U b') = {a Ub) ■ (^6^ U 5® j = (5^ U 6^) ■ {a' U b'). 

As in the construction of the classical Seifert form in knot theory, 

/ (a U 6, a' U b') = I (a, a) + / (a, b') + I (6, a) + / (6, b') : = 

a ■ 6^ + a ■ Sf + b ■ + b ■ Sf = 0 + [a] ■ [b] + [b] ■ [6] = 0 + (a, Qb) + 0 

where the hrst term is zero since the perturbation is towards A, the [ ] denote 
classes in Hk{M), and the “dot” on the last line denotes intersection in M. The 
last term is zero since f4 = 0 G 

Let V' = 6^ U U n consists of a collection of arcs counted by 

a ■ b =2 {a,Qb). Now form V by further perturbing V to be transverse to V in 
X. Each arc oiV AV' becomes a transverse intersection point oiV AV. Thus we 
conclude x ■ x =2 V HV =2 (a, Qb), provided all other points oiVAV come in 
pairs. 

Before we perturbed V toV,Vn V consists of the just described arcs (counted 
by a ■ 6) and a closed 1-manifold S (ZV (Z AU BUC. We denote A U i? U C by T. 
Any additional intersection points arise as the zero locus of a normal perturbation 
to S within the {d — l)-manifold sheet of T in which it locally lives. 



Figure 6.4. 
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But the total number of points along S where the normal held vanishes (as 
it reverses direction) must be even. If it were odd, S would reverse orientation on 
the normal bundle to its sheet in T. But this normal bundle is a direct sum 
of two normal bundles: S within V and S within V'. S must either preserve 
or reverse the orientation of the normal R^ in these two bundles together. So the 
sum will certainly have its orientation preserved. Thus, 


0 =2 X ■ X = 


1/ni/ 


=2 {a, Qb) + I {g} I =2 (a, Qb) + even 


(a,Q6). 


The hrst congruence is from the vanishing of Vk+i G 


□ 


Remark on hypothesis of Lemma 6^. As initially introduced. Figure [6T] portra yed 
two loops I and a on the 2-sphere S"^. But let us now reinterpret Figure 6^ 
as two {S^ X CP^)s, I and a, within S'^ x CP^ by simply taking the Cartesian 
product of the entire diagram with CP^. Thus M is no longer four points but 
{4 points} xCP^. X = F^xCP^ has a nontrivial, low dimensional, Stiefel-Whitney 
class: W 2 {t{X)) is the generator of {S‘^ x CP^) corresponding to (* x CP^) 


under the Kiinneth decomposition. So Lemma |6.4 does not apply, and indeed 


the sign formula in Theorem |6.1| is not valid in this example. Since Wi produces 
isotopic “before 
-|-1; however. 


and “after” submanifolds x CP^, we expect the sign to be 


□ 

Returning to the proof of Theorem 6.1, consider the Mayer-Vietoris sequences 
for P U C, the initial multi-balloon; A U C, the hnal multi-balloon; and AU B, the 
Wilson balloon. We use the abbreviations a fl c, bOc, and a fl 6 for the dimensions 
of the intersected kernels A n C, P fl C, and A n P, respectively. The sequences 
read in part: 


( 1 ) 

P ^ ^ Pfc(AnC) ^Hk{A)®Hk{C) -^Hk{AuC) ^ Hk-i{AnC) , 


( 2 ) 

0 ^ ^ Hk{B nC) ^ Pfe(P) © Hk{C) ^ Pfc(P uC) ^ Pfc-i(P nC) ^ , 

( 3 ) 

0 ^ ^ Hk{A n P) ^ Pfc(A) © Pfc(P) ^ Hk{A U P) ^ Pfc-i(A n P) ^ • 

A basic fact is that the alternating sum of dimensions in any exact sequence with 
held coefficients is zero. We are only interested in collecting mod 2 information so 
we consider just the non-alternating sum of dimensions, let A = s(A U C) + s(P U 
C), the mod 2 “difference” of semi-characteristics. For notational convenience, we 
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extend the definition of the semi-characteristic to odd dimensional manifolds with 
boundary such as A, B, and C, writing s(yl) = 2, where hi is 

again the Z 2 -Betti number. From sequences [T] and we see: 

k k k 

5(^4 U C) = a n c -I- hi{5) + hi{A) -|- hi{C), and 


i=0 

k 


i=0 

k 


i=0 

k 


s{BUC) = h^c + Y, H5) + HB) + h{C), 


so 


1=0 


i=0 


(6.7) 


A = anc-|-6Uc-|-s(A) + s{B), and by Lemma [O 
But from sequence 


i=0 


bk{S) 


cb b 5 ( 2 !) -|- 


( 6 . 8 ) 

But 

(6.9) 


s(A U i?) = a n 6 -f s(A) -|- s{B) + bi{M). 


2 = 1 


+ ^ mod 2 


2=1 


SO 


( 6 . 10 ) 


2 2 

combining lines ( |6.7[ ), ( |6.8[ ), and (6.9) we get 

X{M) 


A = s{AU B) + 


Thus the change A in the semi-characteristic of a fluctuating hypersurface {A U 
C) upon NOT on the support {A U B) of the Wilson balloon operator is given by 
the right-hand side of (6.10), a local formula in agreement with our sign rule. 

□ 


Note 6.5. The computation of x ■ x above relies on a certain regularity of the 
cellular k-chain V representing the class x. We have implicitly used that V is a 
pseudo-manifold, meaning that it is a union of (k + 1 )-cells glued together so that 
every k-face is on the boundary of precisely two (k -\- l)-cells. 


7. Excitations 

Given a GDS on some manifold, and given a ground state T, the state OT may 
be an excited state for certain choices of the operator O. We now briefly discuss 
these excitations. 

One interesting choice of the operator O to consider is a product of Pauli Z over 
all {d— l)-cells that intersect some arc on the dual lattice; we call these dual open 
Wilson loops. This case is completely analogous to the case of the GTG. This 
operator commutes with all terms of the Hamiltonian except those terms near the 
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two endpoints of the arc. The state OT is the same for any two arcs that are 
homologous. The excitations created near the endpoints of the arc must occur 
in pairs: given such a ground state T and such an operator O, there is no state 
whose reduced density matrix agrees with T everywhere except near one of the 
endpoints of the arc, while agreeing with the reduced density matrix of OT near 
that endpoint. 

We can also consider operators obtained from “open balloons”. Assume d = 
2k + 2, k > 0. Choose L to be an “open balloon” (a codimension = 1 submanifold 
with boundary). Let a be the [d — l)-dimensional cellular submanifold supporting 
a basis vector. In this case, choose the operator to be ±NOT, where the sign is 
given by 

Note that in the DS model for d = 2, the operators that drag semions b and 
that drag semions b are related by complex conjugation; however, for the present 
theory, if we replace i by —i in the dehnition of the open balloon operator, then 
the resulting operator is unchanged because n a)) is even. 

Balloon or open balloon operators will either commute or anti-commute with 
dual Wilson loop or dual open Wilson loop operators depending upon the parity 
of the intersection. Further, dual open Wilson loop operators commute with each 
other. The commutation relations between open balloons are more complicated, 
and we leave this as an open question. We also leave the statistics of the excitations 
created by the open balloons as an open question. 


8. Comparison with Z2-Dijkgraaf-Witten theories 


One aspect of the CDS theory (in d = even dimensions) is that topological 
features of the d-manifold on which X it lives may create “missing sectors,” ele¬ 
ments a G Hd-i (X; Z 2 ) so that for hypersurface E C X, [E] = a, 'ip[E] = 0 for 
any zgswf. An example of this is Since y = 2A; -|- 1 = odd and 

Hik-i (CP2^; Z 2 ) = 0, the Hilbert space CDS (CP^^) = 0. 

It is natural to wonder if there is a relation between CDS and the twisted 
Z 2 -Dijkgraaf-Witten (DW) theory, and we thank Meng Cheng and Kevin Walker 
for raising this question and discussions on the resolution. For G = Z 2 , BG = 
K {Z 2 , 1) = RP°°, the inhnite real projective space. 


W {RP°°;R/Z) = 



j even 
j odd 


so there is a unique twisting to consider when j is odd. Generically the effect on the 
Hilbert space of a twisting is to introduce some relations into the corresponding 
untwisted space, so we should check if these relations could cause the “missing 
sectors” observed in GDS. Theorem 18.11 below shows this is not the case. 
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In the untwisted DW theory for a hnite group G, the Hilbert space DW (X'^) = 
is C-linear formal combinations of isomorphism classes of G-principal bundles over 
X. A twist is given by the choice of a cocycle a, [a] G (BG; R/Z). (The class 
[a] determines the theory.) Now, following [10], TDW {X’^) is generated by maps 
/ : A —)■ BG with a relation for every homotopy F : X x I ^ BG, F |xxo= /o and 
F |xxi= fi of the form /o = ufi, where oj = F*{a)[X x /]. From this dehnition, 
it follows easily that if for all maps g : X x ^ BG, g*{a) [X x = 1, then 
dimTDW(A) = dimDW(A), i.e., there are no additional relations. We use this 
criteria to prove: 

Theorem 8.1. For d even, GDS is distinct from Z 2 TDW (and also Z 2 DW, the 
untwisted theory). 

Proof. For simplicity, begin with d = Ak, where GDS is distinguished from both 
twisted and untwisted Dijkgraaf-Witten theory by 

GDS (CP2^) = 0 

and 

Z 2 TDW ^ Z 2 DW (CP^^) ^ C. 

The hrst equality results from y (CP^*^), the Euler characteristic, being odd which 
implies that the empty hypersurface may be “swept over” CP^^ and back to itself 
with an odd number of Morse transitions, each of which counts for (—1) in the 
wave function. 

Since tti (CP^^) = 0, there is only the trivial Z 2 -principal bundle so 

DW (CP2^) ^ C. 

We will prove that any map g : CP^^ x —)■ PP“ homotopically factors through 

the second factor: 


g : CP2^ X 




RP° 


(commutes up to homotopy) 


This immediately implies (by a Z 2 -Stokes theorem) that 

g*{a) [CP^x^^] =g'* =1, 

the second equality, since has no higher homology with any coefficient system. 
Now we construct g'. Fix 6 & and consider 
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The lift I exists since tti = 1. I extends to a map over cone (CP^^) and 

this cone further projects back into RP°°. 


r 

RP^ 

Homotopically, this factors g through gi : S (CP^^) V —)■ PP°°, S denoting 
suspension. But E (CP^^) is also simply connected so gi extends over a cone on 


cone 




PP~ 


Figure 8.1. 


S (CP^^), to a map g 2 , 

g 2 : S {CP^^) y S^U cone (S (CP^^)) ^ PP“. 

But the source of g 2 is homotopy equivalent to 5*^ and making this identihcation, 
g 2 becomes the desired g'. 

With this warmup we can now handle the case d = 4k + 2, k > 1. This is 
done by considering the Hilbert spaces for RP^ x CP^^. The argument will be an 
iterative version of the previous. Since y (PP^ x CP^^) = 2k + 1 = odd, it again 
will suffice to show that for g : RP^ x CP^^ x —)■ PP°° homotopically factors 
through a map g' : RP^ x ^ RP°° (since {RP^ x R/S^) ^ 0). 

To simplify notation, set M = RP^ x S^. Initially we have a (product) bundle: 

Cp2k -^ Qp2k ^ ^ ^poo 

■r 

i-' 

where the homotopy extension h exists since tti (CP^^) = 1 . 

Make h transverse to RP°°~^ C PP“ with Y = h~^ (PP°°“^), and write M = 
Y U P, P = M \ Y . By the previous method, g factors through gi, 

gi : [uy X CP2^ U (5(z/y) x cone (CP^^))] U P ^ PP“, 
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where uy denotes the normal bundle to Y in M. Notice the space in brackets is a 
bundle E over Y with a simply connected hber: S (CP^^): 

S (CP2^)-^ E PP~ 

hi / 

F 

As before, tti (SCP^^) = 1 implies a homotopy factoring through hi. Now make 
hi transverse to RP°^~^ C RP°° and set hi^ (PP°°“^) = L cY. Set Y \ L = Q, 

F = Lug. 

Now Qi factors through g 2 , 

92 '■ \_{^L X S U i9 (z/i) X cone (s (cp^'^))] u g u p ^ pp“. 

Again the space in brackets is a bundle with simply connected fiber: S (S (CP^^)): 

S2 (CP2^)-^ P PP“ 

/ 

/ 

M Y h2 
/ 

L 

Finally, make ^2 transverse to RP°°~^ C RP°° with h^^ (PP°°“^) = J and L\.J = 
K, L = J U K. Now g 2 factors through g^, 

gs : [{uj X (CP2^)) U d (z/j) x cone (S^ (CP^'^))] U P U g U P PP“. 

Again the space in brackets is a bundle with simply connected hber (CP^): 

E3 (CP2^)- G RP^ 

J 

But now J is merely a hnite number of points indexing the components of G. Each 
such component can now be lifted to S°°, coned in S°°, and then projected back 
to RP°°. The result is to factor g^ through g^ where the source of g^ is the source 
of gs with these cones attached to the components of G. But with these cones the 
source of g 4 is homotopy equivalent to RP^ x = J U K U Q U P. Using this 
homotopy equivalence, g^ may be identihed with the map we have been seeking, 
gh □ 

9. Discussion and Relation to Other TQFTs 

We have constructed a Hamiltonian which generalizes the double semion TQFT 
to higher dimensions and we have computed the dimension of the space of zero 
energy states for several different choices of ambient space. For odd dimensions 
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d, we have shown (Remark 5.7) that the ground state subspace of the theory is 
related that of the toric code by a local unitary transformation. 

For even dimensions d, we conjecture that: for all non-empty choices of ambient 
space, there is no local unitary transformation that maps the ground state sub¬ 
space of the GDS theory into the ground state subspace of the GTC or twisted Z 2 
Dijkgraaf-Witten model. More strongly, we conjecture that: given a ball B con¬ 
tained in the ambient space, no ground state wavefunction of the GDS theory can 
be mapped by a local unitary transformation in such a way that the reduced den¬ 
sity matrix on B agrees with the reduced density matrix on B of either a ground 
state wavefunction of the GTC or a ground state wavefunction of the twisted Z 2 
Dijkgraaf-Witten model. To make this conjecture more precise, one may consider 
a sequence of cellulations of the ambient space with bounded local geometry and 
increasing number of cells so that the diameter of the cells tends to zero; then one 
constructs a sequence of quantum circuits which map some ground state wavefunc¬ 
tion of the GDS theory so that the reduced density matrix on B has the desired 
property. These circuits are chosen to minimize the depth of the circuit (taking 
the gates in the circuit to have bounded range); then the conjecture is that this 
depth diverges. 

This conjecture is one way of formally stating that the Hamiltonians correspond 
to different phases of matter. The main evidence for this belief in even dimensions 
is that we have shown that for certain choices of ambient space the ground state 
degeneracy of the GDS model does not agree with that of either the toric code or 
the twisted Z 2 Dijkgraaf-Witten model; this immediately implies that there is no 
local unitary transformation that maps the ground subspace of one theory onto 
that of the other theory in those cases as the dimensions are different. However, 
more work is needed to establish the stronger conjecture above. This is perhaps 
reminiscent of the case in two dimensions where it was shown hrst that no ground 
state of the toric code on a torus can be mapped by a local quantum circuit to 
a product statefH], then it was shown that the ground state of the toric code on 
a sphere could not be mapped to a product state na, and hnally more general 
invariants were found showing that the topological S-matrix of the theory is in¬ 
variant under local unitary transformations [T3] . Hence, to show such a conjecture 
here may require a more detailed analysis of excitations of the theory. This may 
require a description in terms of higher categories |14j. 
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Appendix A. Generic Cellulations 

Consider any set S of points in general position in d-dimensions. Here “general 
position” means that no more than d+1 points he on any hypersphere and that no 
more than n + 1 points lie on any n-dimensional hyperplane for n < d. Consider 
the Voronoi cells dehned by this set of points. Consider any subset T of the d- 
dimensional Voronoi cells. Let B be the boundary of this subset. We now show 
that B is a. d — 1-dimensional PL manifold. 

Let X be any point in B. Then, x is at some distance D from some set of points 
Pi,P 2 , ...,Pk in which are centers of cells in T and points qi, which are centers 
of cells not in S. Thus, there are k +1 points in S at distance D from x. Because 
the points in S are in general position, k + I < d + 1. Further, k > 0 and I > 0 
(otherwise, x is not on the boundary), so k, I < d. 

For notational simplicity, translate the point x to the origin so \pi\ = \qj\ = D. 
The points in B near the origin are the points that are the same distance from 
centers of cells in T as from centers of cells not in T. For any point y, the distance 
from y to a point pi is equal to — 2pi ■ y + y"^. So, a point y near the origin 
(i.e., sufficiently close to the origin that pi, ...,Pk, <?i, ■■■,qi are closer to y than any 
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other point in S) is in B if and only if 

(A.l) maxipi ■ y = max^-gj ■ y. 


Since the points are in general position, the vectors pi, qj span at least k + l — 1 
dimensions. Suppose that the vecto rs Pi ,qi span k + I — 1 dimensions. Then, 
the set of points satisfying the Eq. (A.l) is times the set of points in 

a k + I — 1 dimensional space (i.e., the space which is the span of Pi^qj) which 
satisfy Eq. (A.l). So, in this case we can reduce the problem of showing that B 
is a PL manifold to the problem of showing that B is a PL manifold in this case 
that k + I = d + 1. 


Suppose instead that the vectors v span k + I dimensions. Then, there is some 
nonzero vector v such that pi ■ v = qj ■ v = constant for all i,j. Then, translating 
the origin by this vector n, we arrive at a new problem in which the vectors pi, qj 
span only k + I — 1 dimensions. This reduces to the above case where Pi, qj span 
k + I — 1 dimensions (this vector n is a vector in the space above). So, 

in general we can assume that k + I = d + 1 and the points span d dimensions. 

Consider the space of vectors v such that 


(A.2) V ■pi = v -pj 

for all i,j and 


(A.3) 


vq,^v qj 


for all i,j. These equations give (fc — !) + (/ — 1) = d — 1 constraints. Hence, there 
is at least a one dimensional space of vectors satisfying Eqs. ( A.2[A.3 ). Let n be a 
nonzero vector satisfying Eqs. ( A.2|A.3 ) with v ■ Pi = Zp for some constant Zp and 
V ■ qj = Zq for some constant Zq. By the assumption that the points are in general 
position, Zp ^ Zq, as otherwise there would be d + 1 points on the hyperplane of 
hxed inner product with v. 

There is a d — 1 dimensional space of points orthogonal to v. Let m be any 
vector in this space. Then, 


(A.4) 


satishes Eq. (A.l) if and only if 


y = m + cv 


(A.5) 


c = 


maXjPi ■ m — maxjqj ■ m 

Zq - Zp 


bo, setting y = m -\ - - - ——v gives a PL homeomorphism from a neigh- 

Zq Zp 

borhood of the origin in R'^~^ to a neighborhood of the origin in B. 

Note that without the assumption of general position, B need not be a manifold. 
As an example consider four points in the plane, at coordinates (±1, ±1). Consider 
the set of Voronoi cells containing the cells corresponding to the points at (+1, —1) 
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and (—1, +1). Then, the boundary of this set is the horizontal and vertical axes, 
which has a singularity at the origin. 



